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In the framework of closed time path thermal field theory we investigate the production rate of 
hard thermal photons from a QCD plasma away from equilibrium. Dynamical screening provides a 
finite rate for chemically non-equilibrated distributions of quarks and gluons just as it does in the 
equilibrium situation. Pinch singularities are shown to be absent in the real photon rate even away 
from equilibrium. 



I. INTRODUCTION 



Ultrarelativistic heavy ion experiments at CERN and RHIC focus on energetic direct photons as promising ob- 
servables of the expected QCD plasma phase. Namely the spectrum of thermal photons could be an observable 
effect directly related to the multiparticle dynamics of the plasma system, and as such has been the subject of major 
theoretical interest. As a result the hard thermal photon spectrum has been obtained consistently in thermal field 
theory implementing dynamical screening mechanisms in the framework of hard-thermal-loop (HTL) resummation 
under the assumption that the plasma be in perfect thermal equilibrium |^,^. This later assumption is to be 
contrasted with parallel investigations of the plasma expansion from which the picture of a hot gluon gas |Tl| ] 

has emerged by now: quarks but also gluons are expected to be essentially undersaturated during important parts of 
the plasma history. The effect of this undersaturation on the photon spectrum has been examined in |l^-p^ on the 
basis of folding scaled distributions on conventional matrix elements. The mass singularity appearing in this process 
has been cut-off by hand introducing a thermal quark mass suitably generalized from the equilibrium one. No attempt 
has however been made to obtain the necessary screening from a consistent treatment within thermal field theory so 
that the validity of this cut-off prescription seems to be 'unclear' at the moment. 

In this note we therefore examine in some detail the generalization of the HTL-resummation program to the case 
of undersaturated distribution functions as applied to the hard photon rate. We do so on the basis of the closed time 
path approach for field theoretic systems out-of-equilibrium. We take into account the appearance of recently found 
pinch singularities Jl6| , pT| which can be shown to provide extra contributions to physical rates in out-of-equilibrium 
situations We however demonstrate that these can effectively be neglected in the present case of real photon 

production. Our derivation is also easily applied to the particular case of equilibrium distributions with non- vanishing 
chemical potential 1 19 2^ . Here the photon rate can be explicitly shown to be free of infrared singularities as obtained 
in pc| ] only by numerical integration of the phase space. 

After giving a short review of the theoretical framework used we start with a close look on the contribution from 
the soft part of phase space in section II A| an d add the contribution from the hard part in section 11 B . Finally the 
result is obtained and discussed in section |IIl[ 



II. THERMAL PHOTON RATE 



We adopt the real time formulation of quantum field theory [piy22[ in order to calculate the photon production 
rate. The Keldysh variant is appropriate for a description of systems away from equilibrium p3| . Both propagators 
and self-energies acquire a 2 x 2 matrix structure in this formalism. The components (12) and (21) of the self-energy 
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matrix are related to the emission and absorption probability of the particle species under consideration p3|-|25| . In 
the present context we treat the photons as external probes, i.e., thermally decoupled from the quark-gluon plasma. 
To lowest order they are produced from annihilation and Compton processes 

q + q-> 9 + 1, q{q) + g ^ q{q) + 1, (i) 

but there are essentially no back-reactions that would absorb photons present in the medium. The rate of photon 
emission can thus be calculated as 

from the trace of the (12)-element of the photon-polarization tensor. 




P' 

FIG. 1. The photon polarization tensor ni2 to first order in Us for real photon production. The momentum labels define 
the notation used in the text. External vertices are to be of type 1 and 2, respectively, internal vertices are to be summed over 
a, 6 = 1,2. 



We calculate II12 to leading order in the strong interaction of quarks and gluons. The relevant diagram in fixed 
order perturbation theory is shown in Fig.^ The relevant (12) and (21) propagator components depend on the 
(off-equilibrium) distributions of quarks and gluons. 

Having in mind locally thcrmalized but undersaturatcd distributions for quarks and gluons we neglect higher than 
second order terms in the cumulant expansion for the off-equilibrium correlations . Separating a macroscopic scale 
X from a fast microscopic scale p in the sense of a Wigner transform and applying a gradient expansion to the former, 
the relevant propagators can be shown to depend on a Wigner distribution n{X,p) formerly in the same way as in 
equilibrium |2^ . For the bosonic case all relevant quantities have been listed in [|l8| . Sticking to the conventions used 
there we use the gluon propagator in Feynman gauge as 

iDlliX,p) = -~g^,27reipo)S{p^)n{X,po), iD%{X,p) = iD]^^{X, -p). (3) 

We choose |7|J|,|ll 

n{X v^)-n{v^]-{ Vi3(bo|) Po > , . 

n(X,po) - Hpo) - I _(i + AgnsdpoD) Po < ' 

with riBdpol) the Bose distribution in equilibrium and fugacities Xg parameterizing the deviation from chemical 
equilibrium for gluons. The dependence on the macroscopic variable X is chosen to reside in the parameters Xg = 
Xg{X) and T = T{X). For our approach to be valid we assume that this macroscopic scale is large even against the 
soft scale (gT)^^ for a small QCD coupling g so that we differ in this respect from the assumptions discussed in [ p7[ . 

Fermionic quantities are obtained by substituting fermion distributions n{X,pQ) ~ n{pQ) obtained from (^) through 
the replacement XqUb — > —XqUp, with npipo) the Fermi-Dirac distribution and Xq parameterizing accordingly the 
deviation from equilibrium for the quark- as well as antiquark-distributions: 

"(^'^°)^"(^°) = ll-VHbo|)) Po<0- 
We make use of the (12)-part of the fermion propagator for massless quarks which reads 
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iSuip) = -^'2'^S{p'^)e{po)h{X,po), iS2i{p) = -iSi2{-p)- 



(6) 



All quantitative results will be given for the choice of distributions Eqs.(||) and (||). They thus fit into the framework 
of the model developed in 1^,^ and unnecessary technical complications are avoided because of the factorized depen- 
dence on fugacities. Our main result, i.e. the screening of the quark mass singularity, is however rather independent 
of this choice. It applies in particular equally well to the case of Jiittner distributions 

"■^('P"')" elPol/T±A ' 

used e.g. in the analysis by Strickland |2^. Here the functional dependence on the fugacity models the introduction 
of a chemical potential. However this parameter takes the same value for fermions and anti-fermions, A = A^, and it is 
also attributed to the gluons, A = Ag. This is to be contrasted with the equilibrium distributions, where the chemical 
potential /i enters with different signs for particles and antiparticles, and vanishes for gluons. 

For real photons the self-energy H12 in (^) receives contributions only from spacelike loop-momenta p^ < 0. The 
fixed order result from the diagram Fig.^ turns out to contain IR singular contributions from small — > 0. For the 
equilibrium case it has been shown that this unphysical behavior can be cured taking HTL-contributions from all 
orders of the quark self-energy consistently into account ■ In the spirit of the HTL-resummation program it is 
important to distinguish regions of hard from those of soft momenta |2£| . We choose to separate these scales along 
the line p^ — — fc^ with /c^ to be chosen on an intermediate scale jsj: 

- Plard -T'> kl ^ gT^ > -plf, ~ g'T\ (8) 

For soft momenta resummation of leading HTL-contributions from one-loop self-energy insertions is accomplished 
by substituting the effective resummed propagator as indicated by a blob in Fig.^. The resulting partial rate is IR 
finite and is shown to exactly match the hard contribution thereby introducing the thermal mass-parameter rriq as an 
effective IR cut-off. 




P' 

FIG. 2. The photon polarization tensor 1X12 for soft momentum p ~ gT. The blob indicates the HTL-resummed quark 
propagator. 

Turning to the off-equilibrium situation now our first task is to determine the appropriate generalization of the 
effective resummed propagator to be used on the soft scale. 



A. soft part 



In the non-equilibrium case the analogous resummation program is complicated by the appearance of extra terms 
already at fixed one-loop order. Being proportional to products of basic retarded and advanced propagators 

Ah(p)Aa(p) = , , . —, — r = , 2\2 I — 2' (9) 

-I- i£e[po) p — i£e{po) (p ) + £ 

these terms are in addition plugged with pinch singularities for p^ — ^ [p^ . The resummed propagator for the off- 
equilibrium situation was discussed in for the scalar case with the objective of providing a self consistent cut-off for 
these pinch contributions. Generalizing to the fermionic case resumming successive (one-loop) self-energy insertions 
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is*,{p) = isl^ + isU~^^ab)isl^ + isU-^i:ab)isU-i'^de)isl^ + • • • (10) 

may be rearranged into the effective propagator 

p - h + lepQ p - h* - tepo 

[(l-n(po))(-»Si2)+n(po)(-^S2i)] ^ ~' . , (11) 



^ — T, + iepQ ^ — S* — iepo 

making use of the identities 

/mE = -^(I]i2-S2i), i?eE = i?eSii, (12) 

where S* denotes complex conjugation. 

In the expression Eg. (pi]) ±iepQ is stated exphcitly in order to keep track of the retarded or advanced nature of 
individual terms for later use, even though the complex thermal self-energy E renders this infinitesimal displacement 
superfluous. With the help of Eq.(|l2|) E can be decomposed as 

E = i?eEii + i(Ei2-E2i). (13) 

Evaluating E using equilibrium field theory and restricting to the leading HTL-contribution the first term in 
( pi] ) reduces to the HTL-resummed effective propagator as used in the equilibrium calculations |5|^. The second 
extra term associated with pinch singularities vanishes in the equilibrium calculation because of the detailed balance 
relation E12 = — e^P°/^E2i. It does so as well, for the case of an equilibrium chemical potential /i applied to the 
quark distribution. In the general non-equilibrium situation a non-vanishing self-energy E provides a self consistent 
cut-off for the pinch singularity in Eq. (|ll| ) , as it was already observed by Altherr jl^ . 

We proceed by evaluating the different terms of Eq.(pl]) in the spirit of the HTL-approximation, i.e., assuming the 
external momentum of the propagator to be soft and focusing only on leading self-energy contributions resulting from 
hard loop momenta. First, the imaginary part of E, Eq. (p^) , is linear in the distributions |^ 

-*E-(p)^-^(Ei2(p)-E2i(p)) 

= -g^Cp J ^{^ + ^)5ik^)e{ko)Siip + k f)e{po + ko)[n{ko) + hipo + fco)], (14) 

with a color-factor Cp — 4/3 for Nc — 3 colors. Focusing on the leading contribution only, dimensional and angle 
integrations in Eq.([l^) decouple, 

-i^-{p)\^^^ = -^^ dEkEk{n{Ek)+n{Ek)) j —mP-K), (15) 



with K — (l,fc) a light-like vector with spatial direction of k. In a similar way also the real part of E in Eq.(|le 
can be analyzed in the HTL-approximation. It is found to differ from the imaginary part Eq. (|l^) only in the angular 
integrand. 

The only difference to the equilibrium situation is thus seen to reside in the modified distributions to be integrated 
over Ek = \k\ in Eq . (|l5|) . The combination of these functions present in Eq.([l5|) is important for the following 
discussion. Here it defines the thermal mass parameter mq generalized from the equilibrium case as: 

<-^Cp EdE{n{E)+n{E))^^-—cA\, + ^], (16) 



where in the second step we perform the integration for the specific choice of distributions, Eqs.(^ and (^), respectively. 
There is a factor of 2 in Eq. (^6|) resulting from equal contributions for positive as well as negative /cq to the integral in 
Eq. (|l4[) . In case of an equilibrium chemical potential 2{n{E) -t- n{E)) is replaced by 2nB{E) + np{E — 11) + np{E + /i) 



in the integrand of Eq.(16). 

We now turn to evaluate the second, additional contribution to the propagator in Eq. ([Ill) in HTL-approximation. 
The coefficient in square brackets can be decomposed as 
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(1 - n(po))Si2 + n(po)S2i = (1 - 2n(po))S- + E+, 
where has already been evaluated in Eq.([l^) and 

= 2 (^12 + S21) 

is demonstrated in the following to be non- leading. It contains non- linear terms in the distribution functions: 



:{i> + ^)e{ko)5{k')e{po + ko)5{{p + kf) 



n{po + ko) - n{ko) + 2n{ka)n{po + fco)], 



(17) 



(18) 



(19) 



and keeping only leading contributions the angular integral decouples as before. The accompanying dimensional 
integral turns out to be suppressed by one power of g with respect to the result (|l^ ) for soft external p, po ~ gT and 
hard loop momentum k T: 



dkk[h{po + k)- n{k-po)]{l + 2n(k)) 



d 

dkk — n{k){l + 2n{k)) - 0{g^T^) 



dk 



(20) 



This causes the contribution of E"*" in the resummed propagator Eq.(pT|) being of order 0(1/T) to be dominated by 
the leading contributions from E~ and i?eEii which are of order O {\/{gT)). Neglecting therefore S+ in Eq.([l7|) and 
substituting the remainder for the square brackets in 8*2 Eq.(pl]) the resummed propagator may be rearranged into 
a form similar to the equilibrium case [|6| : 



(70 - ■yp)Im-^ + (70 + jp)Im^ 



D. 



(21) 



Here I)^{po,p), D^{po,p) determining the quark dispersion relations are as in equilibrium pO|-|3^, however, with the 
thermal mass now depending on the off-equilibrium distributions as given in Eq.([l6|). 

It is worth to point out that in equilibrium iS*2(j)) is proportional to np(po), i.e. n{pQ) in Eq.(^) is replaced by 
npipo), which in the soft limit po — is approximated by 1/2. This way the multiplicative factor becomes the same 
as in the non-equilibrium case as it is derived in Eq. (pl|) . 

In order to obtain the partial thermal photon rate from the soft part of phase space we come back to Eq.(^ and 
evaluate II12 from Fig.0 as 



d^p 



Tr [7^ iSl2{p)\HTL 7m*^2i(p -q) + ^tSuiph^ ^S^P - <D\ 



HTLi 



(22) 



with the resummed propagator Eq.(|2l|) and Cq the quark charge. We now follow analogous steps as in equilibrium 
|^,||. Restricting the validity of the present resummed approach by kc ~ ^/gT as discussed above, the factor 1/2 
in Eq.(21) amounts to the approximation npip) — 1/2 applied in the equilibrium calculation. Performing the phase 
space integration for large 3> T we find the soft contribution to the thermal photon rate 



dR 



soft 



2m2(Aq, Ag 



(23) 



As in the equilibrium calculation the modified thermal mass provides a self-consistent IR cut-off under the logarithm. 
It is crucial to note the combination of fugacities present in Eq.(p3|). With the appropriate changes in the thermal 
mass according to Eq.(|l|) the same result is obtained with Jiittner distributions Eq.(Q). For the case of an equilibrium 
chemical potential /i < T ^ E-^ the multiplicative fugacity factor Ag in (^) resulting from the hard fermion line would 
be absent . The partial rate from the hard part of phase space to which we turn next has to match this structure 
in order that all singularities are screened. 
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B. hard part 



For a discussion of the partial rate resulting from the region of hard exchanged/loop momentum it is convenient to 
introduce Mandelstam invariants from the momentum variables defined in Fig.|l| as 

s={k + qf, t^p^ for q + q^g + J, (24) 

and accordingly for the Compton process. In terms of these the kinematic region that remains to be considered is 
expressed as 

-s + kl<pl,rd = t<-kl 2fc2<s<oo. (25) 

Here the calculation is to be done in fixed order e^g^, i.e., from the diagram shown in Fig.^ We need the quark 
propagation with one self-energy correction. For real photon production this amounts to keeping only the second 
term in Eq. ( [lO| ) . The corresponding propagator is obtained by direct calculation of 

iSSM ^ ^St{~^J:ab)^S",2, (26) 



using the relations (13), or alternatively from an expansion of the resummed result Eq. ([ll|) . At 0{g'^) the correction 
to the propagator is found to be 

SSuip) = -HpoKAI - A\)^ReEip)^ - ^fiipo^Al + A\)^{E,2 - ^2i)i> 

-ArAaM'^ - ^(po))Si2 + n(po)S2i]7^. (27) 

Here the first term, which is proportional to 

{Al-Al)^2me{p^)5'{p'), (28) 

corresponds to the virtual correction to the quark propagator and does not contribute to real photon production, 
because the delta function has no support within the relevant region of phase space (^5|) . 

The second term would provide the result found in JpUq], when evaluated in equilibrium. Being proportional to 
ImS it corresponds to a vertical cut of the diagram Fig||7 Switching to Mandelstam variables s and t = the factor 

A| + Ai = 2pl (29) 

appearing here induces a pole in the t-channcl of both the Compton and annihilation process, which is to be regularized 
in the sense of a principal value. Together with one extra factor of t from the spin structure of the quark propagators 
this term is at the origin of the logarithmic mass singularity cut-off by fc^, 

t . fs-P 

dt- 



The third term in Eq.(p7|) contains the pinch singular contribution Eq.(|^). It can be observed though that within 
the region of phase space under consideration, Eq.(p5[), which does not include the zeros of either s, t, u the common 
pole position is approached from t < only. In this situation both pinch and principal value turn out to be equivalent 
as the regularization parameter e can be taken to zero: 



lim / dt- = P / dt— = / dt-. (31) 

As a result associated terms containing the distribution n{po) cancel each other. The one-loop correction to the quark 
propagator therefore reduces to 

SSMp)\,^<-,^^=-j^^^i2A (32) 

which agrees with the corresponding expression in thermal equilibrium up to the fact that E12 is now to be evaluated 
with the more general propagators for the present situation. Calculating 
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j4 

- ffli2 (9) = y j^Tr [7^z5i2 [p - g)7^*<J52i (p) + 5 ^ 6S] , (33) 

we can again rely to a large extent on the work done before in the equilibrium calculation. 

of Eq.(^3|) contains contributions from both annihilation and Compton processes. The phase space integration 
for both partial rates may be written ||33|| 

f°° /■"''^ dt 

/ ds -i^(s,t)$(s,t,^^;A„Ag), (34) 

where the matrix elements F{s,t) are the same as in equilibrium. In <f>(s, t, i?^; Ag, Ag) for the kinematic variables 
a number of integrations can be performed resulting in boundaries Bi,2 for the remaining two dimensional 
integrations: 

-'i>{s,t,EY,Xg,Xg)= [ dEim{Ei) ( dE2n2{E2){l±n3{Ei+ E2- E^)){aEl + bE2 + c)-^^'^ (35) 
s Jb, Jbo. 



s 



with a,b,c functions of Ei. Here the actual combination of distribution functions, rii — n{ni — n), differs for the two 
processes and the sign in the enhancement (blocking) factor appears correspondingly. The expression Eq.(p5|) has 
been evaluated in previous work ||^,|33t assuming Boltzmann approximation for the incoming particles with energies 
Ei,E2- 

In order to obtain an IR finite emission rate it is crucial to keep full quantum statistics for the incoming particles as 
well in contrast to the treatment in |p"^- [l5| , [l9| , p8| . To demonstrate this we first concentrate on the singular contribution 
alone which we isolate with the help of a partial integration over di — —dt/s in Eq.(^). The contribution from the 
t-pole in the annihilation process can be written in the limit kc — > 0: 

It-pole = j ds jln $(s, t = 0, E^- AO - <f (s, t = -s, E^- X,) 



/■°° d 
/ ds di[\nt--t] —.<^{s,t,E^;X, 
Jo Jo dt 



)• (36) 



The logarithmic singularity present in the first contribution is seen to be related to processes with zero momentum 
transfer i = as expected. This restricts the energy of one incoming particle to coincide with the energy of the 
unobserved final particle: 



00 



^ann{s,t^O,E-,) ^ —n{E-,) dEn{E) {1 + n{E)) . (37) 

Js/4E^ 

The same expression is found for the phase space <I>a„„(s,u — 0,Ey) weighting the singular contribution from the 
M-pole in the annihilation process. The Compton process on the other hand contributes two singular terms from quark 
as well as anti-quark scattering each weighted by 

<^>cornis,t^O,E^)^^,,,r^is,u = 0,E^)^^fl{E^) [ dEn{E){l - n{E)) . (38) 

'^1 Js/iE-, 

Summing finally all four singular contributions, products of distribution functions under the remaining integral cancel 
between the two processes leaving behind the same linear combination 2[n{E) + n{E)\ as it appeared in the definition 
of the thermal mass Eq.(|l^). This is crucial for the latter to provide the necessary cut-off. Treating quarks and anti- 
quarks differently for the case of an equilibrium chemical potential the fermionic distributions in Eqs.(^) and ( |3^ ) 
have to be adjusted correspondingly summing again to 2nB{E) + np{E — fj.) + npiE + p) defining the thermal mass 
in this case. On the other hand from Eqs.(|3^) and ( ^ ) it is also evident that applying the Boltzmann approximation 
to the incoming particles spoils this nice feature in either case, and therefore the singularity remains under this 
approximation. 

Interchanging the remaining two integrations over s and Ei and again assuming the photon energy to be large 
with respect to T, i.e., h(Ej) ~ Xqe~^'^/^ , the singular part of the rate in the off-equilibrium situation is now easily 
obtained as 



E 



d^q 



2 3q! 2 -E /T^ ( ^E-fT 



= ^^4^\m^^-"^/Mn^ . (39) 

hard \ '^c 
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In order to complete the calculation we treat the remaining two contributions in Eq. ( |3q ) in order to fix the constant 
under the logarithm in Eq.(|39|). This requires the unrestricted kinematic phase space ^(s, t, E^; Xq, Xg) to be calculated 
with full quantum distributions. This can be accomplished by expanding Bose and Fermi distributions as 



(40) 



For the equilibrium case the resulting extra contributions have been shown to cancel between the annihilation and 
Compton channel |^ making the Boltzmann approximation for the final result better than it is to be expected on the 
basis of individual processes This cancellation, however, depends on the properties of the equilibrium Fermi-Dirac 
and Bose-Einstein distributions and is spoiled for the off-equilibrium situation. For the specific choice Eqs.(U) and 
the additional contributions can nevertheless be evaluated without having to resort to classical approximations. 
We choose to order the resulting terms according to the powers of fugacities present, having in mind that eventually 
Aq <C Ag. Taking out one common factor Xq we obtain for the remaining contribution 



, 2aas 



el^e-^-'/'^T^Xq C{E^,T,Xq,Xq), 



hard 



(41) 



with 



C{E^, T, Xq, Xg) — X 



+A„ 



2 + (l-y)7 + (l-y)lnf -C. 



+ XqXg 



^-y + 4-2)(7 + ln^) + ^C'(2) + ^ln2+(C+-C-) 



(42) 



The appearing symbols are Euler's constant 7, the derivative of the Riemann ('-function evaluated at C'(2) and the 
sums 



^ ^ (_\n 



(43) 



The result Eq.(^) is successfully checked for consistency by noting that the constants present in the equilibrium 
result |6j are reproduced for A, = A^ = 1. 

Summarizing this subsection we find the partial rate for photon emission from the hard part of phase space, Eq. (p5|) , 
to be given by the sum of the two contributions Eq.(^ and (^, 



E 



dR 



_ ^dR'"'3 

hard 



E- 



dR" 



hard 



d^q 



(44) 



hard 



III. RESULT AND CONCLUSIONS 



Calculating the complete emission rate for real photons as the sum of the two partial rates Eqs.(|23D and (|4J) we see 
that the dependence on the arbitrarily chosen fc^ cancels. In the process the generalized thermal mass is established 
as self consistent cut-off for the logarithmic singularity thereby justifying the approach taken in previous work. Our 
result can be written 



dR_ ^aa^ 2 -E,/T 



^(A. + ^)ln 



2E^T 



-\ :tC{E^,T, Xq, Xg 

TT 



(45) 



The dynamical screening of the mass singularity actually does not depend on the explicit form, Eqs.(Q) and of 
the distribution functions assumed in Eq. (p5[) . In fact it only relies on the crucial combination of distributions used 
to define the thermal mass in Eq.(^6|) to appear in front of the singular factor from the sum of Eq.(^7|) and (|3^). 
In particular Jiittner distributions Eq.(|^) lead to a result differing in the thermal mass parameter and the constant 



8 



factors under the logarithm, but the photon spectrum is equally free of singularities. The same conclusion holds also 
for an equilibrium situation with a non- vanishing chemical potential |20| . 

The second important result of our analysis is that we show explicitly the absence of any additional contributions 
from pinch singular terms to the off-equilibrium photon production rate (^5|). As far as the soft momentum scale 
contribution is concerned these terms are shown to be subleading with respect to the dominant HTL-contributions. 
For the hard scale on the other hand the absence of pinch singular contributions is due to the restricted kinematics of 
real photon production and does not apply in general, e.g., for virtual photon production the relevant region of phase 
space is sufficiently enlarged for these terms to become relevant |lj,|22l . 




Turning finally to more quantitative conclusions we note that the rate Eq.(p5|) differs in details from the results of 
both Refs. and but indeed justifies the cut-off prescription employed therein. The discrepancy with respect 
to can be attributed to the fact that distributions for the incoming particles have been treated in Boltzmann 
approximation. We find however that it is necessary to keep track of full quantum statistics in order for the screening 
to be complete. Following we finally obtain a handy estimate by restricting to the idealized situation of a fully 
saturated gluon distribution Ag = 1 with a quark admixture below its equilibrium value by a factor Ag = C = 1/4 g. 
In this situation the photon rate is approximated by 



There is an overall factor 2 with respect to [|12| and a rather large constant c ~ 5 under the logarithm depending on 
the actual value of E^jT . 

According to (|4^) the photon emission rate is reduced by a factor of about 3 /8 with respect to the equilibrium case; 
this reduction is however easily compensated by a higher plasma temperature. 
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^ Under this condition the photon rate given in |u^- [L5| , p8t is dominated by the hard Compton process with the cut-off set by 





(46) 
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